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On the General Solution of the Problem of Disturbed Elliptic 
Motion. By E. Neison, Esq. 

There are two distinct methods of dealing with the problem 
of disturbed elliptic motion. One method is to throw the effect of 
the perturbations directly on the three coordinates, or to suppose 
that the disturbed body moves in an irregular curve which is 

N 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at New York University on May 9, 2015 







1879MNRAS..39..149N 



Mr. Neison, On the Problem xxxix. 3, 


lapproximately of tlie form of an ellipse. The other method is to 
[throw the perturbations directly on the elements of the elliptic 
;orbit, by assuming the body to move in an ellipse which is per¬ 
petually changing in form and position. The first method may 
[be called that of perturbed coordinates, and has the advantage 
[of requiring only three variables,'but in turn the disadvantage of 
“'employing more complex differential equations. The second 
method, which may be termed that of perturbed elements, has 
the disadvantage of requiring six or seven variables, but requires 
a very much simpler form of differential equations. 

What is required for astronomical purposes is the polar 
coordinates of the disturbed body. By making these polar 
coordinates the three variables of the first method, their values 
at any instant is at once obtained. But by the second method 
what is obtained is only the value of the elements of the 
ellipse in which the disturbed body is momentarily moving. 
These must be transformed into an expression for the required 
polar coordinates, and this necessitates extra work. Conse¬ 
quently there is this additional disadvantage in the method of 
perturbed elements. 

Despite the disadvantages which have been pointed out, the 
method of perturbed elements is the simplest and most ex¬ 
peditious way of determining the perturbations of the order of 
the first power of the disturbing forces. For this reason it is 
admirably adapted for employment in the planetary theories, 
where powers of the disturbing forces are rarely required. In 
the lunar theory, however, it is absolutely essential to consider 
the square of the disturbing forces. But this renders most 
complex the known methods of determining the perturbations of 
the elements and necessitates an enormous increase in the labour. 
For instead of there being only six powers and products of the 
three variable coordinates, the second method involves the for¬ 
mation of no less than twenty-eight powers and products of the 
seven variable elements. The calculation of these powers and 
products is laborious and difficult. Therefore, under these con¬ 
ditions, the second method involves many times as much work 
as the first method. 

For these reasons the earlier workers on the lunar theory 
—Damoiseau, Plana, Lubbock, Pontecoulant—employed the 
method of perturbed coordinates. The polar coordinates were 
taken for variables, and the differential equations were simpli¬ 
fied by omitting all quantities of a higher order than it was in¬ 
tended to retain in the approximation, whilst by special artifices 
the amount of work was lessened. Poisson strongly advocated 
the method of perturbed elements, and, at his instigation, Ponte¬ 
coulant employed it in certain parts of his lunar theory, but 
pointed out that to calculate the lunar inequalities by this method 
would be a work of immense labour. Hansen made use of a 
method which is really a modification of the method of perturbed 
coordinates, though he did not throw the perturbations directly 
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on the coordinates. In other respects his system of using his 
Snethod was the same as that employed by his predecessors. In 
d^ll these cases, therefore, special means were taken to obtain a 
Evolution which should be complete to a given order, but which 
liould not be extended to any further order, and depended 
Entirely on the special features of the problem. Lately, Hill, by 
^making use of rectangular coordinates and a modified form of 
the differential equations, has successfully determined the value 
of the perturbations which are independent of the eccentricities. 
His method possesses some generality of form, and is not limited 
to the particular problem. It consists of an expression connect¬ 
ing together the different indeterminate coefficients. It is, how¬ 
ever, an expression which gives a very complex form to the 
complete coefficient of the nth term, and is restricted to the co¬ 
efficient of a small class of easy terms. Moreover, it yields the 
value of the rectangular coordinates and requires extra work 
to transform it into an expression for the polar coordinates. It 
is, however, to be remarked that for the definite purpose which 
its author had in view these are matters of indifference. 

Delaunay, alone, has successfully made use of a method of 
throwing the perturbations on the elements of the elliptic orbit. 
He does not make use of the ordinary method of the variations 
of elliptic elements—that is, the method which has been employed 
by Pontecoulant and Le Merrier in the planetary theories, or 
which was advocated by Poisson for use in the lunar theory. 
Delaunay devised a most ingenious method, which really consists 
in determining the value of the elements which will make the 
disturbing function assume a given form ; a form which enables 
him to directly integrate the resulting differential equations. 
It is a most elegant and ingenious method, but is far more com¬ 
plex than any other, and involves an enormous amount of work. 
It is no exaggeration to say that it involves thrice the amount 
of work which would have been necessary to obtain the same 
result by the use of Pontecoulant’s method. It is, moreover, 
purely a special method, a method depending on a long series of 
particular transformations, and it would be impracticable to 
obtain by its means a general expression for the perturbations of 
the coordinates or elements. On the other hand, Delaunay’s 
method has an important advantage over any method involving 
the use of indeterminate coefficients, as Hill’s or my own ; it is 
possible to separately determine the value of each single term in 
the perturbations without its being necessary to determine any 
other. 

It has long been my opinion that, by choosing suitable 
variables, and by transforming the differential equations in a 
proper manner, it would be possible to obtain a complete general 
solution of the problem of disturbed elliptic motion; that is to 
say, a solution which would be so far general that it would be 
possible to write down the perturbations to the nth order of the 
disturbing forces, and be independent of any particular values of 
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I :the elements. Such a solution would be obtained if it were pos¬ 
sible to form a simple expression for the complete perturbations 
jc^of the order of the ?ith power of the disturbing forces. It seemed 
;|to me that if such an expression could be formed, it would be an 
;!jmportant contribution to the theory of disturbed elliptic motion, 
Snd in particular to the lunar theory. For several years, there- 
" fore, it has been a problem which has engaged much of my 
attention. 

It was not long before it was seen that though it would be. 
possible to obtain an expression for the perturbations of the order 
of the nth power of the disturbing forces, yet it was evident that 
it would be very difficult to prevent this expression becoming 
very lengthy and complex. This difficulty has now been over¬ 
come, and a very simple expression has been found for the per¬ 
turbations of the order of the nth power of the disturbing forces. 

In the present paper I propose to give a short account of the 
principal results, sufficiently full to enable them to be under¬ 
stood ; but the details and minor developments are deferred to a 
subsequent communication, when I propose to apply the results 
to some important points in the lunar theory. 

Let alt be the disturbing function when the elements of the 
orbit have their instantaneous values x, y, z, X, Y, Z, the coefficient 
a being supposed to be a constant, and let the general term of R 
be denoted by 

K cos w = (A) x F y G z H cos { (A) +/X + gY + hZ], (1 ) 

where (A) and (X) are supposed to be independent of the 
elements of the orbit, and to be absolutely constant. 

Suppose JR to denote the value of R when the elements 
have their constant mean values. Then JR will be the same 
function of the mean elements, x, y, 0, X, Y, Z, as R is of the 
instantaneous elements, x, y, z, X, Y, Z, so that the general term 
of R can be -written in the form 

R cos w = (A) x^y G zR cos {(A) + fX+gY + hZ }. (2) 

Let each element of the orbit be considered as the sum of 
its constant mean, value, and either of a single periodical func¬ 
tion which represents its perturbations, or of n periodical func¬ 
tions which collectively represent its perturbations, so that, for 
instance, 

x = x + 8x = x + 5,# + +.+ = x + ( 3 ) 

Then R is to be developed in terms of R and these variables, 
which represent the perturbations of the six elements. 

Let v be taken as the general representative of the three 
elements which have been denoted by x, y, z, and I for the repre¬ 
sentative of the corresponding coefficients F, Gr, H; and similarly 
suppose Y to represent the elements which have been denoted by 
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X, Y, Z and i to represent the corresponding coefficients/, g , h. 
Then it may be supposed that 


v = v + dv = v + 25 m v, V = F+2d m V, ( 4 ) 


where or, W represent the periodical perturbations of the 
elements whose mean values are represented by v and V. Assume 
further that 


‘W“ 1 <5y ~^c v cos y p , 

, -m I -r.ru 

v~ l d m v == 2j p B p cos iv p , 

Wp 


SF = 2C P sin y p , 

F= 3 Jp 4 ~ sin w p , 


( 5 ) 


where m denotes the number of accents to be attached to j 
and J. 

Next by successive direct differentiation 


V s 


d*B_ 

dv s 


r(l) 

r (i-O 


B cos w, 


_ 

dV s ~ 


Si s { cos 


s —.B cos w —sin 5 - ,B sin w 
2 2 


( 6 ) 


Put 


D v = t; 


Dr* 


8F 


d_ 
dV* 


then a known form of expansion gives the expression 


R = B + 8 B ={eZ(X>v+v V )}B 

= {3(1 + 0,+ I [D„? + ... )(i+D f + 2 [D k ] 2 + ... )... }fi. (7) 

From the value assigned to D the following expressions can be 
derived:— 

[D v \ n B cos w = (~) {I (I- 1 ) • . . (I-»+ i)}[oi + o 2 + o 3 +. . . + c n ] n 

cos (w + q' 7i + q"y 2 +.. . + q !i y n ), 

[D F ]«i^co 3 ^» = (i)V 2 "•..^>•}(i)’ , [C, + C 2 + C 3 + .. . + C„]»* 

cos (w + q'y x + q"y 2 +. .. + q n y n ), (8) 

[D v ] m [D r ] Wi ~ n B COS W = jgmgfW+l , . , I (I—i).,,(I— 3 »+l)J 

+ c 2 +. . . + [C m + C m+1 +, . . + 0 TO ] M eos (w + <fy x + q n y 2 +. . . + q n y n ). 


Then putting for brevity 

0k)» = {I (I-i) • • • (I-»»+i)}[o 1 + c a + .. . + c w ]», 
(qiC) m = (*)»faY' . . . q m ] [Cj + C ? +. . . + C m ]» 
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the complete value of tlie general term of B will be 

E cos w = R cos w + - R (I c x + g'tfC,) cos (w + q'y^ +. . . 

... 2 y 

1 /i\ w 

+ I» \2 / ^ ^ + 2®°)" cos + V 2 m 7 *») + • • • (10) 

In tlie preceding expressions the coefficients g- q rf . . . (f 1 only 
take the values plus oi* minus unity, and have been introduced 
for the purpose of clearly showing the changes of sign produced 
by changes of sign in the argument. 

Let 

R = X {B p + Ap]cos w p , (i I) 

where denotes the sum of the coefficients of the argument 
w p which arises from the portion of It depending on the pertur¬ 
bations, so that B p -f A p denotes the complete value of the co¬ 
efficient of the argument w p in the value of It. The function A 
is therefore a constant, and a function of B p , the coefficients I, i, q } 
and the indeterminate coefficients c and C. 

Let (R) be a function of similar form to R, and such that 


(E) = xCF)y(G) Z (H) Ej ( I2 ) 


so that the general term of ‘(R) may be written in the form 

(E) cos w = (A) x F +( r >yG+(G) z H+(H) cos{(A.)+/X + gY + hZ\ , (13) 

and is evidently merely a particular case of the general term (1). 

Its expansion will be identical in form with that of (10). 
Suppose therefore, 

(R) = 2^ (F) y (G) ^ (H) {B p + (A)p} cos w p , (14) 

where (A) p only differs from A p in F, Gr, H being replaced by 
L+(F), G+(G), H+fR) respectively; so that (A) p is merely 
a particular case of A p . 

Suppose the coefficients of the perturbations in the six 
elements of the orbit to be represented by the 67 c indeterminate 
coefficients c 3 . . . c 6h then the function A p can be written in the 
form 

« 20 (r) <? s + 20' (r) (c s c Si ) + 20" (r) (c s c s c s J + . . . + 20^ (r) (e § ... c Sn ) +. . . , (15) 

where 0 (r),. .. <j> n (r ) are known constants, functions of B r , I r , i 
and c s , ... c 8n are particular cases of the series of indeterminate 
coefficients which represent the perturbations of the elements. 
The summation extends to all the combinations of the suffixes r 
and s, ... s n which will yield a term with the argument w p . 

By means of (10) the value of the coefficient A p as developed 
in (15) can be written down to the ^th order, It is merely 
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Inecessary to know the different combinations of the arguments 
w r , yi, ... which will render w r + qy 1 + . . . + cf l y m equal to 
jthe value w p . Suppose there are l such independent combina¬ 
tions of w r , yi, . . . y k , then by giving to s, . . . s rn and r these l 
[series of values, the complete value of (r) (c 9 . . , c Sjn ) as it 
iexists in A p will be obtained. It only remains therefore to 
'determine these l combinations of the arguments. For this 
purpose let 

It = 2 * B m cos w m . 


By direct successive multiplication form the functions 

{R}\ {i?} 3 , . . . {R}n 

neglecting terms of a higher order than necessary for an 
approximation to the given order. Then the different combina¬ 
tions of m in the value of the coefficient of the argument w p in 
the function {B} m will be the l combinations of r, gj, . . . s n in the 
value of %(p m ( r ) (c s . . . c Sm ) as it exists in the coefficient A p . 

In this way the value of A p can be completely determined to 
the nth order, or the 2nbh order as usually computed, so that R 
may be regarded as completely determined to that order. 

Suppose the system of elements represented by v, V to be a 
canonical system, so that the differential equations which give 
the perturbations in the elements may be written in the form 

dsf rfiRi , «t-v • (?V 1 r *c.\ 

-77 = — = —sin?, - 77 - = -y- = v 1 21R cos w. (16) 

dtdV dt dv N ' 

Let it be assumed that 

/ cos . v 7 . 1 sin , . 

sin (»> dt = +- cos («), 


and integrate the system of equations (16). They become 

tr ’ 8t> = 2- - fi? + 8i?}cos w, 8F = 2--{i? + (8i?)sinro}. (17) 
^ to u tv 

Therefore 


^ eos B ^ + ir vw p &1 ‘ }OOS W 

0, sin Wi> = {Ij “ + x r “ ( A U sin w p ■ 


(18) 


In this manner, amongst the 6 k coefficients c 1? . . . c e7 ., there will 
be obtained 6 k independent equations of the form 

°p = rr <*>' W «. + ••• + r W («.... O +. • • }, (19) 

where, for brevity, I is included under i and C under c. 
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This system of 6 h equations amongst the 6 h indeterminate 
[coefficients can be completely solved without difficulty and the 
[values of the 6 h coefficients c u . .. c 6Jc can be determined to any 
[required order. 

; It would not be difficult to transform the preceding equations 
; into others of a still more advantageous form, but this is un- 
" 'necessary at present. In practice it would be found essential to 
have seven series of unknown coefficients, the extra one being 
required for the perturbations of the mean motion. Therefore, 
although the preceding system of equations is very simple, and 
can be written down without difficulty to any required order, it 
still retains the serious disadvantage of requiring seven variables. 
When indeterminate coefficients are employed, this multiplicity 
of variables is the great difficulty of the method of perturbed 
elements, for it involves the solution of a system of jk equations 
amongst 7& unknown quantities, whereas, when the coordinates 
are made the direct variables, as in my “ General Method of 
Treating the Lunar Theory,” it is only necessary to solve a 
system of 3& equations amongst 3 h unknown quantities, and this 
requires far less labour. Moreover, the preceding method re¬ 
quires the entire system of perturbations to be attacked at 
once, and it is not practicable to determine the value of any one 
coefficient by itself. It remains therefore to remove these dis¬ 
advantages by further developing tbe method. 

Instead of supposing It to be developed in terms of the 
variables represented by vv and IV, let it be developed in terms 
of the n times as many variables represented by c\ 7 j, . . . S n v and 
c)i V, . . . 0 n F, and denoted generally by B m v and c m V. As the 
functions represented by and l m V are perfectly analogous to 
“those represented by Bv and cV, the form of the expansion will 
not be altered, and it will merely be necessary to make some 
slight changes in the notation. Put, then, 


=si, (1,-0... (i,-* + .) 


(20) 




and let 2 denote a summation with respect to m, so that, for 
example, 


[s ST (V?) £]*- {I,(I,- 0 • • • a-*+1)} 

([/A+. • •+« -7. ■ - + DW-- •+«] If (*0 
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Then the general term of R can be written in the form 

COS W = B r COS W r + i- ["^Bf (I r jf + g'v J™) -—I COS («J r + g'tt'sj) + • . . 

2 L 1 ™sxJ 

+ - (-)T 2i?r cur+?vC) -Ty- cos ( Wr +2> s )+... ( 22 ) 

Substitute this value with that of (R) in the differential equa¬ 
tions, and on integration they become 

- 1 5v = - i r — B r cos w r + . . . + - v-— x 

v w r v w r + %lqw s 

d + “i y eos («V + 2 " 2 W*) +. . . 

( 2 3 ) 

b F= - L. — i? r sin + . . . + -l r - - — X 

'■ r - r v f w r +^qw s 


•V W r 



( X )JT) ~\ Flr sin ("’»•+ 2 ?2“V> + • • • 


where, as before, the transformation of R into (R) is effected by 
the addition of the quantity (I), where (I) is an absolutely con¬ 
stant quantity, which is the same for all values of r and s. 

Assume that the perturbations of the elements of the order 
of the mth power of the disturbing forces are represented by the 
expressions 

v~ l d m v = j m i R m cos w, V -• J m — B m sin w, ( 24 ) 

then it only remains to determine the value of the coefficients 
f . . . j n and J' . . . J n and the quantities B' . . . A 11 , and the 
right hand side of the equations (23) will become an explicit 
function of known quantities, so that the complete value of the 
perturbations of the elements will have been determined to the 
wth order of the disturbing forces. 

Neglecting the square of the disturbing forces, the equations 
(22) can be written in the form 


v-*d lV = ± i p ±-B p cost«p,. 

V Wp 


s i^=f sin®,, (25) 


where the suffix r has been replaced by p. 


Comparing with (24), 




Bp = B p . 


(26) 


Thus all the terms depending on /, J ; , B' may be regarded as 
known. 
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As far as these terms are concerned, the complete value of 
the term in the disturbing function which has the argument 
w k may be written in the form 

E cos w = Etc cos wjc-t- ^ B Si (I r \ + ~ Rr cos (w r + q'w s ^) + . . . 

+ rj* \ l fis (Ms + sWs) y ~ 1 Rr cos (wv + 'tqws) + ... (27) 

0 w s -J 

with the single condition that the suffixes r, s l9 s 2 , . • . s n which 
have been used to distinguish quantities belonging to different 
arguments are to be given in turn all the different combinations 
of values which will make w r + q'w s +. . . + q n w Sn equal to 
w k . These values will be identical with the l different com¬ 
binations of the suffixes m in the value of the coefficient of the 
term with the argument w t in the function [B} n or w-th power 
of JR. 

The coefficients j", . . . j n and J ", . . . J n can be determined in 
exactly the same manner, and will give rise to similar expression, 
the principal difference being that the simple coefficients i 89 I s will 
be replaced by more complex functions. For the present it is 
unnecessary to write down these terms seriatim, although it can 
be at once done from equations (23). The general expression 
for the terms of the perturbations of the elements of the wth 
order of the disturbing forces will be 


0 1 5 n v — r 

Wfn 


|^i • • • [ b m 

{ g*. ( i r f s +jw £f\.. g ( jrj? +awr) £f m } 


cos (w r +^qw s ) 


KV= SJ 'r—B r 



(28) 


{ g B s (l r j's + gU. + PI /.) J'. • • g S s {l r fs + fX + [I] J?) y 


sin (w r + \qw s \ 


with the conditions that the constants are to take in succession 
all values which satisfy the equations 


W r + (/w s ^ + q'w s ^ +. .. + q b Ws b = w p 

b 1 + b 2 + b 3 +. . . + b m — b ( 2 9) 

b x + 2,b 2 + 3^ +... + mb m — n— I. 

These expressions may seem to be complicated, but this is in¬ 
separable from the terms being given for the nth. order of the dis- 
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turbing forces. As far as the fifth power of the disturbing forces or 
rn u of the ordinary notation, they may be written down in an 
extremely simple form, and even to the tenth order, or m 21 , they 
are neither complicated nor lengthy. The different combinations 
which can be taken by the suffixes r, s 19 . . . s n will be identical 
with the combinations of the suffixes m in the value of the co¬ 
efficient of the term with the argument w p , in the function {R} b or 
6th power of B. 

In this manner, therefore, the complete value of the perturba¬ 
tions of the elements of the orbit may be determined to the nth. 
order of the disturbing forces, and this forms a complete solution 
of the problem of disturbed elliptic motion. The expressions 
which have been obtained for the perturbations admit of being 
transformed into still more advantageous forms. It is necessary, 
however, to defer these extensions and developments to a place 
in a second and more elaborate communication. It may be 
mentioned, however, that by these transformations the com¬ 
plete value of the perturbations of the elements may be made to 
assume the following form :— 


- (mf Sp+ {mi ra . 

+ {ml ([HI, + • • • + PHI/} (SrSs, . . Ss n ) + • • •) cos 


w„ 


5F= ([Q]f^+{[Q]I[Q']i;}i? r 5 Si + ... (■ 

+ {[air (WE,+• • •+[QXf} • • -bo+■ • •)™ 

In this expression R p , R r , B a ■. . . R sn are the coefficients of 
the arguments w pr w r) w g . . . w s in the ordinary disturbing 
function. The functions denoted by [P]J and [Q]J ai are known 
functions of the coefficients J r , I Sm , i r , i Sm . . . and the integrating 
factors w Sm . . . They may be transformed into explicit func¬ 
tions of the form 

>{'((lA„ + 2Vl 8m + (I) i .J ~ ~ } = [Q]I m , 

e s m j 

only under this form the terms depending on high powers of the 
disturbing forces are rendered more complex and lengthy. 

No restriction has been placed on the number of variables 
represented by v and V. In the expression for the general 
term of B it has been assumed that they will be only six in 
number, or x, y , 2, A, Y, Z. These are the number required by 
the lunar theory. The method can be applied without alteration 
to the planetary theory, when there will be twelve variables, 
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19 , y , y\ 2, z f , X, A 7 , Y, Y, Z , Z', the accented quantities re¬ 
ferring to the disturbing planet. It cau be applied, however, to 
■t$ny number of such variables, and absolutely without alteration. 

In the preceding investigation, instead of employing the 
;S}rdinary elliptic elements, use has been made of the six undefined 
Elements whose mean values are x, y, 0, X , Y, Z , and no restric¬ 
tion has been placed on them beyond the single condition that 
they are to be such as will yield a canonical system of differential 
equations. The perturbations of these elements have been com¬ 
pletely determined in explicit series of functions of the form 

p r = (A)Zyv 1 > 

w a = (\)+fsX+g g Y+h 8 Z, (31) 

(Ir*« + ffW ~ = (Prfs + q^sfr) ~ + {&r9s + S^r) ~ + (HA + q&sK) 

The value of the disturbing function has been already deter¬ 
mined by me to the ninth order in terms of the usual elements of 
the elliptic orbit. Two courses can be adopted to render this 
result available for the purpose of the present method ; either it 
can be transformed into an expression in terms of the six 
elements, x, y, z, X , Y, Z , or, vice versa , the expression (31) can 
be transformed into terms of the usual elements. 

Suppose 

Y+Y+Z = L + A + f> + constant = w ri 

where L, A, B denote the mean longitude of the Moon and the 
position of the perigee and ascending node of the lunar orbit. 
Then, as long as this condition is maintained, the arguments of 
the terms in the disturbing function will remain unaltered, with 
the exception that the coefficients /„ g r , h r will be transformed 
into other particular values f k , g k , h k . 

Similarly the coefficient of the term with the argument w r in 
the ordinary disturbing function may be written in the form 

b l a Pk e a]c ‘q I \ 1 + b 2 ar + h s e 2 + b ^ 2 +.. .), 

where F*, G fc , H fc are numerical constants, and assuming 
x = 0, (a, e, y\ y = tf> 2 {a, e, n), z = 0 3 {a, e, y), 

the preceding may be transformed into the value 

vV'A 1+ c ^ 2 + c *y 2 +c ^ z + ■••); 

the only difference being that F r , G r , H r will be different integers 
to F fc , G*, H*.. 

It is evident, therefore, that by properly defining x, y, z, 
either the usual disturbing function can be transformed into 
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! explicit functions of X , Y, Y, a?, y, z, or the equations (31) trans- 
■^fbrmed into explicit functions a, e, g, L, A, B, and their co¬ 
efficients of F fe G fo H*,/*, ^ 

;g; It may be remarked that the preceding method has been 
Already employed to deduce some remarkably interesting and 
[^Important theorems with regard to the terms of long period in 
"the mean longitude. From their form it is obvious that any 
single term or set of terms can be determined by themselves, 
so that the method is peculiarly well adapted for determining the 
terms of long period which converge slowly and require the 
fourth and fifth powers of the disturbing forces to be taken into 
consideration. In a subsequent communication I hope to be 
able to deduce by its aid some important theorems with regard 
to the secular inequalities in the motion of the planets depending 
on the third and fifth powers of the disturbing forces. 

London, 

December 9, 1878. 


Note on the Presence of Particles of Meteoric Dust in the Atmosphere. 

By A. 0 . Ranyard, Esq. 

In a Paper* read before the British Association in 1852, 
Professor Andrews announced that he had discovered particles of 
native iron in the basalt of the Giant’s Causeway. Having 
reduced portions of the rock in a porcelain mortar to a tolerably 
fine powder, magnetic particles were collected by passing a 
magnet several times through the powder. The particles ad¬ 
hering to the magnet were then placed under the microscope 
and moistened with an acid solution of sulphate of copper. On 
some of them copper was deposited in a manner which indicated 
the presence of native iron. Professor Andrews appears to have 
suggested that the particles of native iron may have been de¬ 
rived from meteors which fell when the basalt was in a plastic 
condition. 

In 1867 Dr. T. L. Phipson published a book entitled Meteors , 
Aerolites, and Falling Stars, in which he statesf that he had 
frequently exposed to the wind a sheet of glass covered with 
some transparent mucilaginous substance in order to catch the 
particles of dust floating in the air. He says: “I have found 
that when a glass covered with pure glycerine is exposed to a 
strong wind late in November, it receives a certain number of 
blach angular particles, some three or four of which may be thus 
collected in the space of a couple of hours. The experiment 
being made far in the country, away from the ‘ smuts’ of a town, 
the black particles show themselves all the same. They are, 
however, not soot or charcoal; they can be dissolved in strong 

* Brit. Assoc. Deports for 1852. Part II. pp. 34, 35. f See pp. 229, 230, 
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